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PROGRAM TITLE: RECOIL II BESSEL FUNCTION OF THE SECOMD KIND, 

ORDER ZERO SUBROUTINE (FLOATING POINT, RELOCATABLE) 



1. PURPOSE 

1.1 To compute the Bessel function of the second kind, order zero of a 
floating point argument, found in the A and R registers. 

2. METHOD 

2.1 For < X < 8 the following formula is used: 

N 
N^(X) « I (Gaam * Log^ |) jjx) + ^ A^, 

n«d 

where 2n 

A « f-lf-^ i2^J-»-- (1 ^. i -J. 1 ^ . 4. i) 
n ^"^^T TnTTCnTJ ^J- "^ 2 * 3 n^ 

-12 
and W is the smallest integer for which .A^ < 10 

2.2 For 8 < X < CD 

1/2 H 

V^^V^ liJ H ^""-^^ [^2n ^^" (^ " f ) - ^2n^l°^«(^ - f > 

n«0 

where 00 o 
iS^... (2^-1)2 

g «* • 3 

^ iiC8xr 

and where N is the first integer for which 



\^X > =2N °^ »2N*1 < ^^""^ ' 



whichever happens first. 
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3. RESTRICTIONS 

3.1 Range of argument X > 

3*2 Contents of L and V loops are destroyed. 

3»3 The following subroutines must be in memory: sin-cos (AN-012.1), log 
(AN.O37), Jq(AN.055). 

3,k The following modifications should be made in AN- 055 so that it will 
use AN-012.1. 



xage 



■a 9 



L (origin + 005^ TRA 



0000 



FCA 
N0P 



7762 



k. USB 

U.l Calling Sequence 



Ot 


1 FCA 
TRA 


ARG 
^0 


a +1 


PZK 
F7,?, 


L( AN-012.1) 
L(AN-055) 


a + 2 


P7.R 
ERR0R 


L(AN-037) 
RETURN 


a + 3 1 


N0HMAL 


RETURN 



h.2 The routine occupies words 0000 thru Q163q or II6 full words. 

Z1-.3 Minimum accuracy: 8 decimal places. When speed is desired and 
several fewer accurate places are needed^ it is suggested that 
10" be entered into location OI62-3, where P is one greater than 
the number of accurate decimal places srequired. 

k.k After a noxmal exit the value of the function is found in the A and R 
registers. Error return is caused by negative arguments 
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5. CODING INFORMATION 
5.1 Constants 

5.1.1 FIoasMng Point 



1.0 


in 


0002 


2.0 


ill 


OOOU 


8.0 


in 


0010 


2/it « 0.63661 9772 


in 


OC^O 


Gainma « 0.57721 566^^9 


in 


0056 


itA «« 0.78539 6l63is- 
10-^^ 


in 
in 


0102 
0162 



5.1.2 Fixed Point 

1 at B 39 in OOO3 
1 at B 38 in OOC^ 

5.2 Ers'.sable locations: 0006-7, 0<05b~5> 0066"7> 0076-7; 0100-1, 0136-7; 

0156-7- 

5.3 ^Sister ISb.'pq includes: 

1) msic roatine, 0000-0163. 

2) Relocation matrix, i{.000-i^l63. 

3) M-OOJ^- Relocation routine, 772^1-7757. 
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O 



Set up 
Exit Addresses 



X : 



< 



4t. 



X : 8 



<8 



>8 



n « 
SUM « 



Compute 



Jfc 



SUM « SUM 4> C 



n[«x] 



azH*! 




X 



> 1 



i 1 




< 1 




n a n-j-l I 

* 




-12 

"xcf I .„.,,„.,....„{^.,i 



3- 



jSUM - I (r + log^ |)Jo^X) 




